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An  earlier  letter  to  the  editor  [Piquette  and  Van  Buren,  J.  Acoust.  Soc.  Am.  79,  179-180 
(1986))  commented  on  an  article  by  D.  Censor  [J.  Acoust.  Soc.  Am.  76,  1527-1534  ( 1984)]. 
In  his  reply  [J.  Acoust.  Soc.  Am.  79,  181-182  (1986)],  Censor  made  some  observations  that 
warrant  further  comment.  The  present  letter  considers  in  detail  plane-wave  scattering  by  a 
vibrating  plane  and  thereby  demonstrates  analytically  that  Censor’s  theory  predicts  sum  and 
difference  frequency  component  waves  that  are  of  the  order  of  magnitude  of  pseudosound.  This 
example  also  illustrates  that  the  parametrically  growing  sum  and  difference  frequency 
component  waves  generated  nonlinearly  overwhelm  those  predicted  by  the  Censor  theory 
within  a  fraction  of  a  wavelength  of  the  scatterer’s  surface. 

PACS  numbers:  42.25.Gf,  43.25.  Jh,  43.25.Cb 


Recently,  we  published  a  letter  to  the  editor1  which 
commented  on  an  article2  published  by  D.  Censor.  In  his 
reply,3  Censor  made  several  points  which  we  believe  require 
further  comment. 

Censor’s  remarks3  seem  to  imply  that  we  are  denying 
the  existence  of  the  Doppler  effect  in  acoustics.  Censor3  has 
pointed  out  that  the  Doppler  effect  has  previously  been  ques¬ 
tioned  by  Petzval.  Petzval’s  objections  were  resolved  by 
Mach  (see  reference  in  Tolman4).  Unfortunately,  Petzval’s 
objections  (and  Mach’s  subsequent  discussion)  appear  to 
bear  little  relation  to  the  questions  we  have  raised1  in  regard 
to  Censor’s  work.2  According  to  Tolman’s  summary,4  Petz¬ 
val  denied  the  existence  of  the  Doppler  effect,  claiming  that 
the  transmitted  frequency  of  a  source  equals  the  received 
frequency  regardless  of  the  physical  properties  of  the  medi¬ 
um,  provided  that  the  continuity  condition  is  satisfied  at  all 
points  at  rest  with  respect  to  the  source,  and  provided  that  all 
particles  of  the  medium  have  identical  velocity  vectors.  We 
agree  with  Mach’s  conclusion  that  Petzval’s  argument  is 
generally  invalid.  Our  arguments  concerning  Censor’s  work 
in  no  way  contradict  the  existence  of  a  frequency  change  due 
to  a  relative  motion  between  source  and  observer. 

We  agree  that  additional  frequency  components  are 
generated  by  the  scattering  of  an  acoustic  wave  by  a  vibrat¬ 
ing  surface  via  the  Doppler  effect.  However,  we  have 
shown5*  that  such  additional  frequency  components  are 
also  generated  by  the  nonlinear  interaction  of  the  primary 
^  waves  in  the  fluid  surrounding  the  scattering  surface  and 
H  that,  in  general,  the  waves  generated  by  the  nonlinear  effect 
overwhelm  those  generated  by  the  Doppler  mechanism 
o  within  a  fraction  of  a  wavelength  of  the  scatterer’s  surface. 

We  wish  to  emphasize  once  again  that  we  do  not  ques- 
UJ  tion  the  validity  of  Censor’s  work  in  the  area  of  electromag- 
*" — ^  netism.  Maxwell’s  equations  are  strictly  linear,  and  nonlin- 
tJL.  ear  effects  enter  electromagnetic  calculations  only  through 
the  constitutive  relations.  For  most  media  the  constitutive 
relations  for  electromagnetism  are  linear  or  essentially  so, 
and  nonlinear  effects  are  negligible.  However,  the  funda¬ 
mental  equations  of  acoustics  that  correspond  to  Maxwell’s 


equations  in  electromagnetism  (i.e.,  the  equations  of  contin¬ 
uity  and  momentum  conservation)  are  intrinsically  nonlin¬ 
ear.  In  addition,  all  fluid  media  are  nonlinear  through  their 
constitutive  relations;  i.e.,  pressure  is  a  nonlinear  function  of 
condensation.  Most  fluids,  such  as  water  and  air,  are  signifi¬ 
cantly  nonlinear.  The  combined  effects  of  these  two  sources 
of  nonlinearity  in  acoustics  can  produce  significant  contri¬ 
butions  to  sum  and  difference  frequency  components 
through  mixing  during  propagation  of  multifrequency 
waves.  Also  in  acoustics,  a  vibrating  boundary  is  a  radiator 
of  sound,  thereby  producing  a  wave  at  the  frequency  of  mo¬ 
tion  of  the  boundary  which  then  interacts  with  the  incident 
and  reflected  waves.  There  is  no  such  radiation  effect  in  the 
electromagnetic  case  (assuming  the  boundary  is  un¬ 
charged). 

We  have  previously  stated1,5  that  Censor’s  theory  pre¬ 
dicts  results  that  are  of  the  same  order  of  magnitude  as  pseu¬ 
dosound.  Pseudosound  is  an  effect  that  arises  due  to  the  un¬ 
certainty  in  the  motion  of  the  acoustic  sensor.  A  sensor  that 
moves  freely  with  the  particles  of  the  medium  measures  in 
the  Lagrangian  frame  of  reference.  A  sensor  that  is  fixed  in 
space  measures  in  the  Eulerian  frame  of  reference.  Acoustic 
problems  solved  in  either  reference  frame  give  equivalent 
results  to  first  order  in  Mach  number.  (Here,  Mach  number 
is  taken  to  be  the  ratio  of  the  maximum  speed  of  the  particles 
of  the  medium  to  the  speed  of  acoustic  waves  of  infinitesimal 
amplitude.)  However,  the  predicted  results  diff'r  to  second 
and  higher  order  in  Mach  number.  Since  it  is  unlikely  that 
one  knows  the  state  of  motion  of  one’s  sensor  sufficiently 
accurately  to  distinguish  to  what  extent  it  moves  freely  with 
the  fluid  (and  thus  measures  in  the  Lagrangian  frame)  or  to 
what  extent  it  remains  fixed  (and  thus  measur-s  in  the  Eu¬ 
lerian  frame),  pseudosecond  must  be  considered  to  be  the 
practical  lower  limit  of  the  measurability  of  acoustic  quanti¬ 
ties.  Censor’s  theory  predicts  sum  and  difference  frequency 
component  waves  that  are  of  the  order  of  magnitude  of  pseu¬ 
dosound. 

Consider  a  simple  quantitative  example:  the  scattering 
of  a  normally  incident  sinusoidal  plane  wave  (of  angular 


frequency  co)  from  an  infinite  plane  vibrating  uniformly  and 
harmonically  (at  angular  frequency  fl).  (This  case  has  been 
considered  in  detail  elsewhere/')  In  order  to  demonstrate 
that  Censor’s  theory  predicts  results  that  are  of  the  order  of 
pseudosound,  we  will  solve  this  problem  in  Lagrangian  co¬ 
ordinates  assuming  that  the  linear  wave  equation  is  exact. 
We  will  then  transform  this  solution  into  Eulerian  coordi¬ 
nates.  Next,  we  will  compare  the  transformed  solution  with 
that  computed  using  Censor’s  theory. 

Censor  has  chosen  to  solve  the  problem  of  the  scattering 
of  an  acoustic  wave  by  a  vibrating  obstacle  in  the  Eulerian 
(spatial)  reference  frame.  The  incident  plane  wave  in  Euler¬ 
ian  coordinates  is  given  by  PE  (x,t)  =  P0eh,,x/C  +  0 .  We  turn 
off  all  nonlinearities  so  that  the  linear  wave  equation  applies. 
The  displacement  of  the  plane  is  described  by  the  function 
£(0 ,t)  —  e  sin(flr).  Applying  Censor's  theory  to  this  situa¬ 
tion,  we  obtain  the  following  expression  for  the  sum  and 
difference  frequency  pressure  components: 


(1) 


PoC-fl 

where  Pr  =  ip0cCle  =  pressure  amplitude  of  the  directly  ra¬ 
diated  wave,  p0  =  fluid  density,  c  =  infinitesimal  amplitude 
wave  speed,  and  0 1  =  CL  ±0. 

Next,  we  consider  the  same  problem  in  Lagrangian  co¬ 
ordinates.  In  this  case  the  incident  wave  is  given  by 
P1  ( a,t )  sb  P0e'a(a/C  f  Since  the  boundary  is  fixed  relative 
to  a  Lagrangian  frame  of  reference,  the  incident  wave  re¬ 
flects  as  from  a  rigid  body  giving  a  reflected  pressure 
P0e~‘“,,a/C~'>.  yhe  racjiation  from  the  moving  boundary 
provides  a  pressure  wave  _ 'n,0/c  ~  These  two  waves 
travel  without  interaction  when  the  nonlinearities  are  turned 
off.  Thus,  no  sum  and  difference  frequency  components  are 
present  in  the  solution  obtained  in  a  Lagrangian  frame  of 
reference.  However,  sum  and  difference  frequency  compo¬ 
nents  arise  if  we  transform  the  solution  into  Eulerian  coordi¬ 
nates.  An  arbitrary  Eulerian  function  fE(x,t)  i;  trans¬ 
formed  into  the  associated  Lagrangian  function  fu  ( a,t ) 
using  the  expansion7 


fL(a,t)  =fB(x,t)  |x  =  a  +  f(fl,f) 

+  df-Ef:x,t)  §{a,n  +  ---, 

OX  x  »  a 

where  n  =  Lagrangian  coordinate,  x  =  Eulerian  coordinate, 
and£  =  displacement  (common  to  both  systems).  Applying 
this  transformation  to  the  Lagrangian  scattered  wave 
P0e  '"J(a/c  0  we  obtain  the  following  expression  (accurate 
to  second-order  Mach  number)  for  the  sum  and  difference 
frequency  pressure  components  that  arise  in  the  Eulerian 
reference  frame: 


P  Wr  - 
P±  = - —00 


*  \{x/c  ~  t) 


(2) 


Note  that  since  P ±  is  zero  in  Lagrangian  coordinates  when 
the  nonlinearities  are  turned  off,  Eq.  (2)  represents  pseudo¬ 
sound.  The  pressure  represented  by  Eq.  (2)  woulu  only  be 
measured  by  an  acoustic  sensor  fixed  in  space  (i.e.,  in  an 
Eulerian  frame),  and  only  if  the  nonlinearities  are  turned  off. 
Since  any  real  acoustic  sensor  must  move  with  the  molecules 


of  the  fluid  to  some  extent,  the  pressure  given  by  Eq.  (2) 
would  never  be  measured  by  such  a  sensor  even  assuming  the 
linear  wave  equation  to  be  exact.  Any  real  acoustic  sensor  in 
such  a  situation  would  measure  a  pressure  between  zero  and 
that  given  by  Eq.  (2),  due  to  the  uncertainty  in  the  hydro¬ 
phone’s  motion. 

Comparing  Eq.  (1)  (the  solution  using  Censor’s  the¬ 
ory  )  to  Eq.  (2)  (which  results  from  merely  transforming  the 
monofrequency  rigid-body  scattered  wave  from  Lagrangian 
to  Eulerian  coordinates  and  retaining  only  the  sum  and  dif¬ 
ference  frequency  components),  we  see  that  they  are  clearly 
of  the  same  order  of  magnitude  (differing  only  in  that  0  now 
appears  in  the  coefficient  in  place  of  0  ±  ).  Thus  the  effect 
predicted  by  Censor’s  theory  arises  essentially  from  the 
transformation  between  coordinates  and  not  from  medium 
nonlinearities.  An  acoustic  sensor  measuring  in  the  Lagran¬ 
gian  frame,  i.e.,  moving  with  the  fluid,  would  not  detect  sum 
and  difference  frequency  components  arising  from  the  effect 
Censor  predicts.  An  acoustic  sensor  measuring  in  the  Euler¬ 
ian  frame,  i.e.,  fixed  in  space,  would  detect  a  signal  due  to 
Censor’s  effect  but  this  signal  would  usually  be  substantially 
smaller  than  that  arising  from  nonlinearities. 

Now  consider  this  problem  from  the  point  of  view  of 
solving  the  nonlinear  wave  equation.  The  sum  and  difference 
frequency  components  arising  from  nonlinearities  for  plane- 
wave  scattering  from  a  vibrating  plane  are  obtained  in  Ref.  6. 
These  can  be  approximated  by 

„  ( ( 1  +  y)P»Pr ^  ,  ,  —  iifij ^ t(x/c— o 

f  *  H"  4W>  '  J  '  <3> 

where  y  is  the  one-dimensional  nonlinearity  parameter, 
1+B/A. 

Note  in  Eq.  (3)  the  presence  of  the  distance  factor  x, 
characteristic  of  parametric  sources,  in  the  coefficient  pre¬ 
ceding  the  complex  exponential  function.  This  distance  fac¬ 
tor  causes  the  nonlinearly  generated  signal  to  increase  in 
amplitude  with  increasing  propagation  distance.  The  energy 
required  for  this  parametric  increase  in  amplitude  of  the  sum 
and  difference  frequency  component  waves  is  provided  by  an 
energy  drain  from  the  primary  waves.  Equation  (3),  ob¬ 
tained  via  a  perturbation  solution,  remains  valid  as  long  as 
the  energy  densities  of  the  secondary  waves  remain  small 
relative  to  the  energy  densities  of  the  primary  waves.  (See 
Ref.  6  for  a  more  complete  discussion  of  this. ) 

It  is  clear  from  comparing  Eqs.  (1)  and  (3)  that  the 
nonlinearly  generated  signal  will  become  greater  than  that 
generated  by  Censor’s  Doppler  mechanism  at  some  propa¬ 
gation  distance.  We  demonstrate  below  that  this  propaga¬ 
tion  distance  is  always  a  fraction  of  the  longest  wavelength 
involved  in  the  problem.  The  distance  factor  x  of  Eq.  (3) 
arises  because  every  element  of  fluid  between  the  scatterer’s 
surface  and  the  observation  point  acts  as  a  source  of  sum  and 
difference  frequency  waves.5  Thus  the  greater  the  propaga¬ 
tion  distance,  the  greater  the  volume  of  fluid  between  the 
scatte.er  and  the  observation  point  and,  generally,  the 
greater  is  the  amplitude  of  the  nonlinearly  generated  signal 
Note  that  the  solution  to  the  nonlinear  wave  equation  does 
not  simply  produce  a  “small  correction"  to  the  solution 
based  on  the  linear  wave  equation,  as  might  be  expected,  but 
rather  predicts  sum  and  difference  frequency  components 
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that  rapidly  grow  to  a  significant  amplitude.  Such  a  possibil¬ 
ity  exists  whenever  two  primary  waves  of  differii  g  frequency 
are  simultaneously  present  in  a  fluid.  The  parametric  growth 
of  the  nonlinearly  generated  sum  and  difference  frequency 
component  waves  generally  occurs  whenever  there  is  a  sig¬ 
nificant  "collimation  region,”  i.e.,  whenever  there  is  a  signif¬ 
icant  volume  of  fluid  through  which  the  wave  vectors  of  the 
primary  waves  point  in  the  same  direction. 

When  we  compare  the  nonlinear  volume  effect  solution 
for  the  sum  and  difference  frequency  components,  given  by 
Eq.  ( 3 ) ,  to  the  solution  for  these  components  resulting  from 
the  boundary  Doppler  effect  predicted,  by  Censor,  given  by 
Eq.  (1),  we  obtain 

±  ]  Censor  surface  4 C _ 

]  nonlinear  volume 

For  a  planar  surface  frequency  of  100  kHz,  in  a  water 
medium,  the  two  effects  become  equal  at  a  distance  *  on  the 
order  of  0. 14  cm.  One  should  not  be  disturbed  by  the  pres¬ 
ence  of  the  ft  term  in  the  denominator  of  the  ratio  in  Eq.  (4) . 
In  the  case  of  low  ft,  the  relevant  factor  to  scale  the  distance 
is  the  wavelength  associated  with  ft,  which  is  2j rc/ft.  Hence, 
if  we  let*  =  d(2~c/fl),  we  can  determine  the  fraction  dot  a 
wavelength  at  which  the  nonlinear  volume  effect  overtakes 
Censor’s  surface  effect.  This  occurs  for*  =  0.09  (2rrc/ft). 
Therefore,  even  in  the  limiting  case  in  which  the  frequency  of 
vibration  of  the  planar  surface  approaches  zero  ( maximizing 
Censor’s  boundary  effect  relative  iO  the  nonli.  r  volume 
effect),  the  sum  and  difference-frequency  pressure  compo¬ 
nents  generated  by  nonlinear  volume  effects  exceed  those 
produced  by  the  boundary  effect  within  a  propagation  dis¬ 
tance  less  than  the  longest  wavelength  involved  in  the  prob¬ 
lem.  Thus,  as  previously  mentioned,  we  do  not  deny  the  exis¬ 
tence  of  the  new  spectral  components  generated  by  the 
Doppler  mechanism  cited  by  Censor.  It  is  simply  that  the 
amplitudes  of  the  frequency  components  generated  by  this 
mechanism  are  generally  much  smaller  than  those  generated 
by  the  nonlinear  mechanism.  Note  that  the  frequency  com¬ 
ponents  created  by  the  Doppler  effect  are  extremely  small  in 
the  present  case  due  to  the  oscillatory  nature  of  the  motion  of 
the  vibrating  boundary.  If  the  boundary  were  instead  mov¬ 
ing  with  constant  velocity,  the  familiar  and  much  stronger 
Doppler-shifted  reflected  wave  would  be  obtained. 

Since  Censor’s  theory  predicts  results  that  are  of  the 
order  of  magnitude  of  pseudosound,  his  theory  would  not  be 
very  useful  in  monitoring  the  motion  of  vibrating  boundaries 
even  if  an  experimental  arrangement  were  devised  that  was 
unfavorable  to  the  parametric  growth  of  sum  and  difference 
frequency  component  waves  generated  by  the  nonlinear 
mechanism.  (Such  a  possibility  arises,  for  example,  if  the 
incident  plane  wave  in  the  problem  discussed  above  arrives 
at  non-normal  incidence. )  As  previously  mentioned,  the  un¬ 
certainty  in  the  motion  of  a  real  acoustic  sensor  is  sufficient 
to  inhibit  the  meaningful  measurement  of  acoustic  quantities 
predicted  to  be  of  the  orde,  of  magnitude  of  pseudosound. 

In  his  reply1  to  our  comments,'  Censor  suggests  using  a 
perforated  scatterer  in  an  attempt  to  create  a  situation  which 
enhances  the  Doppler  mechanism  over  the  nonlinear  mecha¬ 
nism  for  generating  sum  and  difference  frequency  pressure 


waves.  We  have  argued' "  that  since  both  theories  depend  on 
the  Mach  number  in  the  same  way,  any  change  in  the  param¬ 
eters  of  the  proolem  would  change  the  predicted  values  of 
each  theory  by  the  same  factor.  (This  fact  was  originally 
noted  by  Rogers.")  Censor  states1  that  such  perforations 
would  eliminate  the  displacement  of  the  surrounding  fluid 
caused  by  the  vibrating  boundary.  However,  if  the  displace¬ 
ment  were  to  entirely  vanish,  the  quantity  Pr  of  Eqs.  ( 1 )  and 
(3)  would  become  zero,  and  thus  both  theories  would  pre¬ 
dict  null  results.  Of  course,  such  perforations  would  not  en¬ 
tirely  eliminate  the  radiated  wave.  While  an  exact  solution 
for  this  problem  would  require  re-solving  both  theories  for 
the  perforated  surface  involved,  it  is  possible  to  obtain  ap¬ 
proximate  solutions  for  each  theory  by  slightly  modifying 
Eqs.  (1)  and  (3).  This  is  possible  if  the  perforation  diame¬ 
ters  and  spacings  are  sufficiently  dissimilar  from  either  pri¬ 
mary  wavelength  (thus  avoiding  linear  diffraction  phenom¬ 
ena).  The  primary  waves  are  essentially  affected  in  two  ways 
by  using  a  perforated  scatterer:  (i)  The  amplitude  of  the 
primary  scattered  wave  of  angular  frequency  a>  and,  (ii)  the 
amplitude  of  the  directly  radiated  primary  wave  of  angular 
frequency  ft,  are  reduced.  Although  the  amplitude  of  the 
incident  plane  wave  of  angular  frequency  a>  is  not  affected, 
this  primary  wave  produces  very  little  contribution  to  the 
sum  and  difference  frequency  component  waves  generated 
by  the  nonlinear  mechanism,  due  to  the  unfavorable  propa¬ 
gation  direction  of  this  wave  relative  to  the  directly  radiated 
wave  of  angular  frequency  ft.  In  fact,  Eq.  ( 3 )  ( which  is  only 
approximate)  omits  the  small  contribution  to  the  nonlinear 
field  due  to  direct  interaction  of  the  incident  wave  with  the 
radiated  wave,  and  only  accounts  for  the  interaction  of  the 
linearly  scattered  wave  with  the  directly  radiated  wave. 
(The  complete  solution  is  given  in  Ref.  6.)  We  can  account 
approximately  for  the  influence  of  the  perforations  on  the 
primary  scattered  wave  by  adjusting  the  quantity  P0  appro¬ 
priately.  Similarly,  we  can  account  approximately  for  the 
influence  of  the  perforations  on  the  directly  radiated  wave  by 
adjusting  the  quantity  Pr  appropriately.  However,  note  that 
the  quantity  P  ±  ,  as  computed  using  the  Censor  theory,  giv¬ 
en  by  Eq.  ( 1 ),  depends  on  the  quantities  P0  and  Pr  in  exactly 
the  same  way  as  the  expression  for  P±  derived  using  nonlin¬ 
ear  theory,  given  by  Eq.  (3).  Thus  both  computed  quantities 
change  by  the  same  factor  if  the  vibrating  plane  is  perforated. 
This  means  that  the  ratio  of  the  computed  values,  given  by 
Eq.  (4),  is  unaffected  by  the  presence  of  perforations  in  the 
scattering  surface. 

We  agree  with  the  final  observation  Censor  made  in  his 
reply1  to  our  original  letter.'  The  experiment  we  performed" 
could  be  improved  if  an  arrangement  were  possible  that  pro¬ 
duced  sufficiently  wide  nulls  in  the  primary  field  to  reduce 
the  degrading  effect  of  hydrophone  self-nonlinearity  on  the 
measurements.  (We  were  unable  to  devise  an  experimental 
arrangement  to  achieve  this.)  However,  as  previously  men¬ 
tioned,  the  meaningful  measurement  of  sum  and  difference 
frequency  w  aves  predicted  to  be  only  of  the  order  of  magni¬ 
tude  of  pseudosound,  as  is  the  case  with  Censor’s  theory, 
would  not  be  practical  even  using  a  completely  linear  hydro¬ 
phone.  Since  the  Censor  theory  predicts  results  that  are  of 
the  order  of  magnitude  of  pseudosound,  as  we  have  demon- 


1535 


J.  Acoust.  Soc.  Am.,  Vol.  80,  No.  5,  Novemoer  1986 


Letters  to  the  Editor 


1535 


1 flE?  WSEfKrararoCTHKreSlC^XSJnSCaGOGTiraOTBSJIHX^OEI^fiaanifiKH  PJ4  JUi  ■,L?IT.-V- l'.*  ^UltlTF.L'VL^VX'VITn.^ 

1 


strated  above  and  elsewhere,5,6  we  conclude  that  his  theory  is 
not  very  useful  for  monitoring  the  motion  of  moving  boun¬ 
daries  in  the  acoustical  case. 

‘J.  C.  Piquette  and  A.  L.  Van  Buren,  "Comments  on  ‘Harmonic  and  tran¬ 
sient  scattering  from  time  varying  obstacles  [J.  Acoust.  Soc.  Am.  76, 
1527-1534  ( 1984) ],’  ”  J.  Acoust.  Soc.  Am.  79, 179-180  ( 1986). 

2D.  Censor,  “Harmonic  and  transient  scattering  from  time  varying  obsta¬ 
cles,”  J.  Acoust.  Soc.  Am.  76,  1527-1534  ( 1984). 

3D.  Censor,  “Reply  to  ‘Comments  on  Harmonic  and  transient  scattering 
from  time  varying  obstacles  (J.  Acoust.  Soc.  Am.  76,  1527-1534 


(1984)],’  ”  J.  Acoust.  Soc.  Am.  79,  181-182  ( 1986). 
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